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ABSTRACT 

We examine eight young stellar objects in the OMC-2 star forming region based 
on observations from the SOFIA/FORCAST early science phase, the Spitzer Space 
Telescope, the Herschel Space Observatory, 2MASS, APEX, and other results in the 
literature. We show the spectral energy distributions of these objects from near-infrared 
to millimeter wavelengths, and compare the SEDs with those of sheet collapse models 
of protostars and circumstellar disks. Four of the objects can be modelled as protostars 
with infalling envelopes, two as young stars surrounded by disks, and the remaining 
two objects have double-peaked SEDs. We model the double-peaked sources as binaries 
containing a young star with a disk and a protostar. The six most luminous sources are 
found in a dense group within a 0.15 x 0.25 pc region; these sources have luminosities 
ranging from 300 L Q to 20 L . The most embedded source (OMC-2 FIR 4) can be fit 
by a class protostar model having a luminosity of ~ 50 L Q and mass infall rate of 
~ 10~ 4 M & yr- 1 . 

Subject headings: stars: formation — infrared: stars — circumstellar matter — radiative 
transfer — protoplanetary disks 
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INTRODUCTION 



OMC-2 is a star forming region embedded in the molecular filament stretching northward 
of the Orion Nebula and is one of the most luminous regions in the Orion A molecular cloud - 



Recent distance estima tes place it at ~ 420 pc (jHirota et alj|2007l . iJeffries 



2007 



Menten et al. 



2007 



Sandstrom et al.ll2007l ). As a site of star formation activity, it contains a large number of protostars 
and pre-main sequence stars with dusty circumstellar disks ([Peterson &: Megeathll2008l ). Therefore, 
OMC-2 represents a laboratory for the investigation of the detailed properties of protostars and 
disks which are forming in an active environment. 



Gatlev et al.l (|1974l ) first i dentified OMC-2 as a c omplex of pre-main sequence objects associ 



ated with nearby CO emission. iThronson et al.l (|1978l ) measured the infrared luminosity of OMC-2 
to be 2100 L , but did not resolve the far-IR components of IRS 3 and IRS 4 and thus could not 
provide a measurement of their individual S EDs. Later work e xpanded the coverage of OMC-2 to 
higher resolution at the longer wavelengths. IChini et al.l (119971 ) conducted a millimeter study that 
resolved several cold (20 K) dusty cores t hroughout OMC -2/3. In OMC-2, these cores are associated 
with submillimeter continuum emission (ILis et al.lll998l ). The detection of centimeter sources wit h 
the VLA established many of the cores as bona fide protostars (jReipurth. Rodriguez fc Chinilll999l ). 
Recent studies have included high resolutio n, ground-based observations to further examine bot h 
the stellar component and the dusty cores (jNielbock. Chini &: Mullerll2003l . IShimaiiri et al.ll2008l ). 
These studies revealed several stellar sources associated with the dusty cores and substructure in 
the cold dust of FIR 4. 

While previous studies include high resolution observations at 10 and 20 /mi, submillimeter, 
and millimeter wavelengths, there is a lack of data at high spatial resolution in the ~ 20 — 200 
//m range which is critical for constraining models of protostellar envelopes and for characterizing 
any external heatin g of both dusty disks and envelopes. During its early science phase, SOFIA 
dYoung et alJboidl observed the Orion Nebul a Cluster and OMC-2 with the FORCAST mid-IR 
camera (j Adams et al.l [2J31CJ, iHerter et al.l I2O120 and detected eight intermediate-luminosity proto- 
stars and disks in OM C-2. We combine these data with data from t he 2MASS near-infra red survey 
(jSkrutskie et al.1 120061 ) , the Spitzer Orion Molecular Cloud Survey (jMegeath et al.1 120051 ) , the Her- 
schel Orion Protostar Survey (HOPS; Fischer et al. 12010 ) . submillimeter data from APEX, and 
other published data in the literature to produce complete SEDs of these objects from near-IR to 
millimeter wavelengths. Finally, we compare the SEDs to those of model protostars and disks in 
order to derive the physical properties of these objects. 
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OBSERVATIONS 



2.1. SOFIA/FORCAST Observations 



SOFIA/FORCAST observed OMC-2 during the SOFIA short science phase on December 1 and 
December 8, 2010, at an altitude of ~ 13100 m. Two pointings were used: {05 h 35 m , -05° 08') and 
(05 h 35 m , -05° 15'). Each pointing covered 3.4' x 3.2', with a (rectified) pixel scale of 0.77" /pixel. 
The first pointing was observed at 19.7 fim and 37.1 /im (with the dichroic beamsplitter). The 
second pointing was observed at 19.7 //m and 31.4 j iva (with the dich roic), and 37.1 //m (without 
the dichroic). Observing modes C2N and C2NC2 (jHerter et al.1 120121 ) were used. A large (4.8') 
asymmetric chop throw was used to chop off nearby nebulosity and minimize off-axis coma in the 
on-source beam. The chop frequency was 2 — 5 Hz and the nod amplitude was approximately 20'. 
A 5-point, 10" step dither pattern was performed to reject bad pixels during post-processing. The 
integration time in each dither position was approximately 30 sec. 



The data were processed with the reduction steps described in iHerter et al.1 ([2012]). Point 
source photometry was performed on sources in the images using two-dimensional Gaussian fits to 
their radial brightness profiles. Typical FWHM image quality values were 3.3" and 4.3" at 19.7 
and 37.1 fim, respectively. During the December 8 flight, some elongation of the PSF was seen in 
the cross-elevation direction, resulting in image quality of 5.9" and 4.2" in the cross-elevation and 
elevation direction, respectively, in the second pointing field at 31.4 /xm. Photomet ric signal-to- 
noise was consistent with the mean point source sensitivity over 3 short science flights (jHerter et al 



20JJ). 



The extracted photometric fluxes were calibrated to a flux density using a standard instrument 
response derived from measurem ents of standard sta rs and solar system objects over several short 
science and basic science flights (jHerter et al.ll2012l ). We applied small (< 2%) color corrections 
to each source, since the standard instrument response represents the response to a flat spectrum 
{yFv = constant) source. The estimated 3 a uncertainty in the calibration due to variations in flat 
field, water vapor burden, and altitude is approximately ±20%. 



2.2. Spitzer Observations 



Spitzer/IRAC images were taken in two epochs in 2004 March 09 and October 10 as part of 
program PID 43. The reported photometry is the median value obtained from four 10.4 second 
high dynamic range mode frames. The magnitudes of the sources were ex tracted with ape rture 
photometry from the Basic Calibrated Data using the IDL aper.pro program (jLandsmanlll993l ). A 2 
pixel radius and 2 to 6 pixel background annulu s was used (Megeath et al. in prep.), and a standard 
calibration was adopted (see Gutermuth et al. 120091 and Reach et al. 120051 ). The MIPS scan map 
was made in 2004 March 20 with the fast rate and 160" cross scan offset as part of program PID 
58, resulting in 30 seconds of integration time per pixel. The data were reduced with the MIPS 
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instrument team's Data Analysis Tool (Go rdon et al. 2005 ), and photometry was obtained using 
a the IDL implementation of DAOPHOT (|Landsmanlll993l ). DAOPHO T was modified to i gnore 



saturated pixels (Kruykova et al. submitted). The calibration of lEngelbracht et al 



(2007 ) was 



adopt ed. The magnit udes were computed to flu x densities using the zero fluxes given in lReach et al 
(|20Q5l ) for IRAC and Engelbracht et all (|2007l ) for MIPS. Given the relatively insensitivity of the 
fluxes to the spectral shape, color corrections were not applied. 

Spit zer/IRS spectra were acquired for two sources (MIR 29 and MIR 31+32; Nielobock et 
al. [2003) in the SOFIA/FORCAST fields-of-view. These 5 — 36 /iin spectra were obtained with 
the short-low and long-low modules (PIDs 30859 and 50734). We follow the method described in 
Poteet et al.l (|201ll ) to extract and calibrate the spectra. 



2.3. Herschel Observations 



The Herschel Space Observatory observed an 8' square field encompassing HOPS 59, 60, 108, 
368, 369, and 370 on 2010 September 28 (observing day 502, observat ion IDs 1342205228 and 
1342205229) in the 70 /im and 160 lira bands of the PACS instrument (jPoglitsch et al.l hold ) as 
part of the HOPS Open Time Key Program. These data have angular resolutions of 5.2" and 
12", respectively. The target field was observed with homogeneous coverage using two orthogonal 
scanning directions and a scan speed of 20"/s, with each scan repeated 8 times for a total observation 
time of 3290 s. To facilitate point source photom etry, a high-pass filter was applied to the data to 
minimize the effects of drift ([Fischer et al.ll2010l ). 



The aperture photometry was obtained using a 9.6" aperture at 70 fim and a 12.8" aperture 
at 160 jim. with a background annulus extending from the aperture limit to twice that value in 
both channels. The results were corrected with the encircled energy fraction provided by the PACS 
consortium (priv. comm.). The error estimate, which is dominated by calibration uncertainties, is 
estimated to be 5% at both wavelengths. 



2.4. APEX Data 

We include 350 and 870 ^m phot ometry acqui r ed at APEX with SABOCA and LABOCA, 



respectively (proposal ID 088.C-0994). IStanke et al.1 ([2010) describe the observing and data reduc- 



tion procedures for the APEX data. The beam sizes of the final reduced maps are 7.3" and 19" for 
the SABOCA and LABOCA observations respectively. We measured the flux density per beam at 
the position of the source initially using the coordinates from the HOPS survey 70 fim point source 
photometry. If the source is well-detected, we re-centered the coordinates and measure the flux 
centered on the source in order to account for possible pointing offsets between data sets. Other- 
wise, we adopt the PACS 70 //m coordinates and measure a flux which we can consider an upper 
limit on the source flux at the wavelength of interest. The errors are dominated by the calibration 
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uncertainty and are therefore adopted to be 40% of the measured flux at both wavelengths. 



3. RESULTS 

Fig. [T] shows a Spitzer /IRAC and SOFIA/FORCAST multi-wavelength false color image of 
the northern FORCAST field. In both FORCAST fields, SOFIA/FORCAST detected most point 
sources (except the faint source HOPS 63, a likely protostar, and HOPS 64, a likely young stellar 
object), that were present in the Herschel/P 'ACS 70 and 160 /jm images, plus the source SOF 6. 
This field shows a relatively compact group of luminous young stellar objects, with seven sources 
found in a 0.13 x 0.25 pc region. The two remaining sources, SOF 7 and SOF 8 are found to the 
south of the group and are not shown. In Table [fl we list the eight SOFIA/FORCAST sources 
detected in the field and cross-correlate this list with previous identifications. In Table [2j we present 
the flux densities for each source. We used these flux densities to construct a SED for each source 
(Fi g. EL with addition al values at 1-2 /m i (2MASS), 10.4 /mi, 11.9 jum, and 17.8 /mi (Nielbock et 



al. |2003j), and 1.3 mm (IChini et al.lll997h . Fig. [2] also shows that the Spitzer/IRS spectra for SOF 



7 and SOF 8 are in excellent agreement with the FORCAST data. 



4. MODELS 



The physical properties of the SOFIA-detected objects were derived using models consisting of 
an infalling envelope and/or a star+disk system. For the envelope, we adopted the sheet collapse 
models developed by Hartmann, Calvet & Boss (|1996l . hereafter HCB96), which seem appropriate 
given the filamentary nature of the OMC-2 molecular emission. In these models, the envelope 
is flattened not only in the inner region, because of rotation, but also at large scales due to the 
natural elongation of the cloud. The inner region of the envelope is also partially excavated by an 
outflo w cavity. We adopt ed a dust opacity law that is based on the observed class I object L1551 
IRS5 (jOsorio et al.1 120031 ) . The dust mixture includes graphites, silicates, troilites and water ice 
compounds with a standard grain-size distribution of n(a) oc a -3 ' 5 , with a minimum and maximun 
size dust grains of 0.005 /mi and 0.3 fim, respective ly. The temperature distribut i ons a nd SEDs 
were calculated using radiative codes developed by IKenyon. Calvet fc Hartmann! (|1993l ) and by 
HCB96. In some cas es, external heating of the envelope was added by using a single temperature 
modified blackbody (jPoteet et al.ll201ll ) of temperature ~ 30 K; this is expected in the OMC-2 



region which is near the Orion Nebula and is likely heated by the massive stars in that regions. 

For the disk component o f the model, we used the flared accretion disk models developed 
by iD'Alessio et al.l (119991 . 120061 ). In these models the disk density is given by the conservation 
of angular momentum flux and an uniform mass accretion rate. The temperature distribution is 
determined from heating by viscous dissipation and stellar irradiation. The disk radius is assumed 
to coincide with the centrifugal radius (which is the largest radius on the equatorial plane that 
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receives the infalling material). The disk models also account for dust evolution given by grain 
growth and settling towards the disk midplane. We assumed that the dust composition and grain- 
size distribution is the same as in the envelope, except in the disk midplane where a maximum 
grain size of 1 mm is adopted. 

The fits were found by v isual inspection, following the SEP beha viour as a function of d i fferen t 
physical parameters given in lDe Buizer. Osorio fc Calvetl (|2005l ) and iD'Alessio et al.l (|1999l . 12006). 
In some cases, the inclination angle could be constrained from available near-IR images. 



5. DISCUSSION OF OBJECTS 

The best-fitting model SEDs are overlayed with the data in Fig. [2j The constrained model 
parameters for seven objects are listed in Table [3j We provide comments on each object: 



SOF 1: The A > 70 /im fluxes of this source are considered upper limits due to contamination from 
extended emission associated with SOF 2N. Therefore, we fit this object by taking into account 
only the SED at A < 70 /xm. The best fitting model consists of a young star surrounded by a disk 
with a high mass accretion rate (5 x 10~ 6 Mq yr _1 ). 



SOF 2N: ICohen fc Frogell (|1977l ) presented near-IR images showing that IRS 4 has two distinct 
components, which are flagged as confused i n the 2MASS catalog. Both components (MIR 21 and 
MIR 22) were resolved by Nielbock et al. (|2003) at 10.4 and 17.8 /xm. We assume the northern 
component emits most of the flux at A > 17.8 /xm, as evidenced by the relative flux densities at 
17.8 /xm, and this is the component we model. A high luminosity of 300 Lq as well as some external 
heating is required by our model to explain the SOFIA flux at 37 /im. Therefore, SOF 2N likely 
hosts an intermediate-mass protostar. A highly flattened, nearly edge-on envelope with a large 
opening-angle cavity is also required to reproduce the short wavelength data. 



SOF 3: This deeply embedded source is a class protostar associated with FIR 4 and was not 
detected by Nielbock et al. (2003) at 10.4 ^m. This object has a large envelope mass (2.5 - 10 
Mq) and and dominates the submillimeter emission in this field. Two fits were obtained: one 
which treated the 160 micron as a detection, and one which treated 160 um as an upper limit 
(which would be the case if the 160 um emission were contaminated by the extended emission). 
The resulting luminosites range 30 — 50 Lq and the implied infall rates r ange from 4 x 10~ 4 to 
1 x 10 -4 M & yr _1 ; the highest infall rates in our sample (see Table 3). IShimajiri et al.1 (2008) 
resolved several small dust clumps toward FIR 4 at 3.3 mm, but none of these clumps shows a peak 
that is clearly coincident with our 8 — 160 /iin detections. 



SOF 4: The SOFIA data reveal that the SED has a double peak. Our attempts to model this 
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source as a single object failed to explain the two peaks. Rather, we interpret this source as a dou- 
ble stellar system, where one of the stars (more evolved) is surrounded only by the disk component 
that dominates the emission at A < 40 /iin. The other star (likely in an earlier evolutionary phase) 
is deeply embedded into its infalling envelope that dominates the emission at A > 40 /im. The best 
fit to the disk source requires a high mass accretion rate (M acc = 10 -5 M Q yr _1 ). 

SOF 5: Bipolar nebulosity from this object is seen at 3 — 8 /im (Megeath et al., in prep.). This can 
be used to constrain the inclination angle, resulting in a better estimate of the luminosity and the 
density inferred from the modeling. 

SOF 6: This is a young stellar object with a disk. It was not detected at A > 30 /xm and we do not 
model this source. 



SOF 7: The model for this object requires a central luminosity of 15 Lq and a contribution from 
external heating. A Spitzer/IRS spectrum shows a deep silicate absorption at 9.7 /jm which is fitted 
by assuming a high inclination angle and a large cavity. Inspection of the 3.6 and 4.5 /mi images 
shows a small lobe of nebulosity extending ~ 4" to the northeast. 



SOF 8: Nielbock et al. (|2003l ) resolved a double source (MIR 31 and 32) with projected separation 
2.4" at this location at 10.4 and 11.9 /im. This object has a nearly flat SED from ~ 2 /im to 160 
/im. We achieved a best fit to the SED using a binary system model. The best fit model indicates 
that in one of the sources the emission arises from a disk around an intermediate mass star and 
produces a luminosity of 130 L & . The companion was fit by a protostar model that contains a 
dense envelope of 3 Mq irradiated by a central luminosity of 7 Lq. 



6. CONCLUSION 

We have used SOFIA early science observations to examine the detailed properties of eight 
intermediate-luminosity (7 — 300 Lq) protostars and circumstellar disks in OMC-2. Based on the 
model fits, it is believed that four of these objects are protostars and two are young stars with 
circumstellar disks, while two are likely young binary systems. This work demonstrates the poten- 
tial of SOFIA to observe protostars in active, luminous star forming regions, while simultaneously 
complementing both space- and ground-based observations. With such observations, we can char- 
acterize in detail the properties of protostars and examine the dependence of star formation on 
environment. 
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Fig. 1. — False color image showing the northern SOFIA/FORCAST pointing field in OMC-2 at 4.5 
fjan (blue, IRAC), 19.7 fim (green, FORCAST), and 37.1 /an (red, FORCAST). SOFIA/FORCAST 
sources 1-6 are labeled; see Tabled] for previous identifications. The protostellar candidate HOPS 
63 was not detected by SOFIA/FORCAST. The remaining background point sources seen at 4.5 
//m are dominated by photospheres and young stars with disks. 
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Fig. 2.— SEDs of SOFIA/FORCAST sources SOF 1 - 8. The data points are color coded by data 
set (facility). Uncertainties are shown with error bars and downward arrows indicate data points 
that are upper limits. Solid black curves: Total SEDs for the best-fit models. Dotted green curves: 
Disk component of binary systems. Dashed red curves: Envelope-dominated component of binary 
systems. Dotted black curve: Lower luminosity model fit for SOF 3. Note that the absorption 
feature at 9.7 /im of source SOF 7 shows some emission, likely because of contamination by nearby 
sources. Only the red side of the feature has been fitted. The SED shortward of 11 fim is dominated 
by noise and has not been included in the fit. 
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Table 1. Positio ns and identification s for all sources detected by SOFIA/F ORCA ST (SOF). IRS 
IDs are given in[Gatlev et al] dl g74|), MIR IDs are giv en in Nielbock et al. j2 Qoj), FIR IDs 



given m 



Chini et al. 



arc 



Lis et al.l (|1998l L 



(jl997h and lShimaiiri et al.l (j20od ). CSO IDs are given m 
and VLA IDs are given in Reipurth et al. (|1999i ). Positions are derived from the Spitzer/IRAC 
3.6 /im image, except SOF 2N (derived at 19 //m) and SOF 3 (derived at 8.0 |im). 



SOF 


HOPS 


RA (2000) 


DEC (2000) 


IRS 


2MASS J 


MIR 


FIR 


CSO 


VLA 


1 


GG 


05:35:26.84 


-05:09:24.7 


2 


05352683-0509244 


20 








2N 


370 


05:35:27.57 


-05:09:33.8 


4N 


05352762-0509337 


21 


3 


22 


11 


3 


108 


05:35:27.08 


-05:10:00.3 








4 


23 


12 


4 


369 


05:35:26.98 


-05:10:17.4 


3 


05352696-0510173 


27 


5 


24 




5 


368 


05:35:24.72 


-05:10:29.9 


1 


05352477-0510296 


28 






13 


6 




05:35:26.89 


-05:11:07.6 




05352686-0511076 










7 


GO 


05:35:23.39 


-05:12:02.6 




05352341-0512023 


29 


6b 


25 




8 


59 


05:35:20.14 


-05:13:15.7 




05352014-0513156 


31+32 


6d 


30 





Table 2. Flux densities and uncertainties in Janskies for all sources detected by SOFIA/FORCAST' 



F4.5 



F:,.. 



Fg.O 



F19.7 



F2-1 



P31.4 



^37.1 



F70 



F ieo 



F350 



F 8 70 



SOF 1 0.532 ± 0.001 0.879 ± 0.002 2.321 ± 0.003 3.059 ± 0.004 4.86 ± 0.97 5.31 ± 0.26 



SOF 2N 0.330 ± 0.002 0.829 ± 0.005 0.831 ± 0.005 
SOF 3 



15.0 ± 3.0 45.24 ± 10.42 
1.519 ±0.031 



0.030 ± 0.001 

SOF 4 1.304 ± 0.002 2.145 ± 0.004 3.958 ± 0.004 5.844 ± 0.006 10.8 ± 2.2 17.51 ± 0.66 

SOF 5 0.330 ± 0.002 0.921 ± 0.002 1.423 ± 0.003 2.486 ± 0.003 19.1 ± 3.8 25.58 ± 0.17 

SOF 6 1.345 ± 0.001 1.476 ± 0.002 1.457 ± 0.001 1.490 ± 0.005 1.24 ± 0.25 1.291 ± 0.010 
SOF 7 0.018 ± 0.001 0.047 ± 0.001 0.044 ± 0.001 0.044 ± 0.003 •■■ 1.79 ± 0.01 

SOF 8 0.968 ± 0.001 2.890 ± 0.004 4.034 ± 0.001 5.040 ± 0.004 10.5 ±2.1 11.2 ± 0.1 



4.48 ± 0.90 < 59.5 ± 3.0 < 202 ± 10 < 12.8 ±5.1 < 2.9 ± 1.2 
227 ± 45 718.1 ± 8.1 482 ± 17 26.7 ± 10.6 < 5.6 ± 2.2 
8.4 ±1.7 121.2 ± 4.9 468 ± 35 < 44.5 ± 17.6 < 7.9 ± 3.2 



13.9 ± 2.8 
65 ± 13 



6.9 ± 1.4 
16.3 ± 3.3 



12.3 ± 2.5 
23.5 ± 4.7 



48.0 ± 3.2 
104.8 ± 4.0 



64.8 ± 1.5 
55.4 ± 1.7 



206 ± 23 
105 ± 9 



89.9 ± 10.1 
69.7 ± 10.8 



20.0 ± 8.0 

5.9 ± 2.4 



10.0 ± 4.0 

4.3 ± 1.7 



< 4.8 ± 1.9 
< 1.41 ± 0.56 



< 2.14 ± 0.86 

< 0.88 ± 0.35 



a Scc also Peterson Mcgcath ( 2008) for Spitzer photometry of sources detected at 24 iim. 
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Table 3. Constrained model parameters for FORCAST-detected OMC-2 sources. a 



Envelope Disk 



Source 


L 


i 


V 


R c 


6 


Renv 


PlAU b 


M„f b ' c 


Monv 


-Rdisk 


Mdisk 




(Le) 


(degrees) 




(AU) 


(degrees) 


(AU) (10- 


13 g cm~ 3 ) 


(M e yr- 1 ) 


(M Q ) 


(AU) 


(M @ ) 


SOF 1 


GO 


40 
















40 


1 


SOF 2N 


300 


80 


3 


300 


40 


5000 


1.5 


~3xl0~ 5 


0.8 


280 


1.6 


SOF 3 


50 


50 


2 


380 


8 


5000 


20 


~4xl0~ 4 


10 


380 


0.6 


SOF 3 d 


30 


70 


2 


380 


8 


5000 


5 


~lxl0~ 4 


2.5 


380 


0.6 


SOF 4 C 


20 


70 


2.5 


100 





5000 


9.0 


~2xl0~ 4 


4 








60 


40 
















40 


1.5 


SOF 5 


40 


■50 


2.5 


100 


■5 


10000 


1.5 


~3xl0~ 5 


2 


100 


0.2 


SOF 7 


15 


70 


2.5 


100 


30 


5000 


2.0 


~4xl0~ 5 


1 


100 


0.8 


SOF 8 C 


7 


■50 


2 


100 





5000 


6.0 


~lxl0~ 4 


3 








130 


30 
















100 


0.1 



a L is the total luminosity of the source, i is the inclination angle of the polar (rotational) axis of the 
system (that is measured with respect to the line of sight), 77 is a measure of the degree of flattening of 
the envelope (1 < r\ < 4, with 77 = 4 corresponding to a very flattened envelope), R c is the centrifugal 
radius, 9 is the aperture angle of the envelope cavity, i? e nv is the outer radius of the envelope, piau is a 
reference density (corresponding to the density at a radius of 1 AU of a spherically symmetric free-falling 
envelope with the same mass infall rate), Mi n f is the mass infall rate (see note c), M cnv is the envelope 
mass (obtained by integration of the density distribution) , i?disk is the outer radius of the disk, and M^^ 
is the disk mass. 

b Values are likely upper limits because of possible contamination by nearby sources of FIR/submm flux 
densities. 



Estimated from eq. 3 of Kenvon. Calvet fc Hartmanu ( 1993h . assuming a nominal value of 1 M© for 



the stellar mass. M; n f is proportional to piAV and to the square root of the stellar mass, which is not 
well-constrained. Values are accurate within a factor of two for stellar masses in the range 0.25-4 Mq. 

d Fit obtained when F\qq is taken as an upper limit. 



c Modeled as a double source. In one of the sources the emission is dominated by a star+envelope and 
in the other source it is dominated by a star+disk. 



